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Inhomogeneity-Induced Superconductivity? 
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A t- J- like model for inhomogeneous superconductivity of cuprate oxides is presented, in which 
local anisotropic magnetic terms are essential. We show that this model predicts pairing, consis- 
tent with experiments, and argue how the macroscopic phase-coherent state gradually grows upon 
lowering of the temperature. We show that appropriate inhomogeneities are essential in order to 
have significant pair binding in the thermodynamic limit. Particularly, local breaking of SU(2) 
spin symmetry is an efficient mechanism for inducing pairing of two holes, as well as explaining 
the magnetic scattering properties. We also discuss the connection of the resulting inhomogeneity- 
induced superconductivity to recent experimental evidence for a linear relation between magnetic 
incommensurability and the superconducting transition temperature, as a function of doping. 



There is a growing body of experimental evidence sug- 
gesting that the superconducting state in cuprate oxides 
is "inhomogeneous," such that the locally defined charge 
density varies across the sample in the ground state. 
Spatially inhomogeneous features in the spin and charge 
channels have been indicated in a number of experiments 
on high-T c materials The simplest realization of 

this state is the so-called "stripe" phase where charges 
cluster in nanoscale linear patterns and the remainder of 
the sample is essentially an antiferromagnetically corre- 
lated insulator. This represents a nanoscale distribution 
of charge and spin, rather than a global phase separation. 
These experiments lead us to a central question: Is the 
superconducting state found in high-T c cuprates inhomo- 
geneous as a result of spin/charge inhomogeneities? We 
believe that the answer to this question is yes. Moreover, 
we argue that spatial spin/charge inhomogeneities are in 
fact necessary for pairing and subsequent formation of 
the superconducting state in these compounds. This sit- 
uation should be contrasted with the case of conventional 
superconductors, resolved by the BCS (Bardeen-Cooper- 
Schrieffer) theory, that starts with a homogeneous metal- 
lic state and describes the formation of a homogeneous 
superconducting state. It is commonly believed that 
magnetic correlations, characterized by the spin exchange 
energy J ~ are responsible for the pairing inter- 

actions in the cuprates and are, therefore, crucial for our 
inhomogeneous exchange approach. Moreover, the exis- 
tence of a spin gap has been experimentally proven 0. 
A model that naturally incorporates these features (inho- 
mogeneities, magnetism and spin gap) is a t- J- like model 
with explicitly broken spatial and magnetic symmetries. 

Of central importance for the present work is a micro- 
scopic model which captures the main low energy physics 
of doped antiferromagnetic (AF) Mott insulators. In par- 
ticular, we show that our minimal model properly de- 
scribes the magnetic properties observed in a wide va- 
riety of doped cuprate oxide materials. The key ingre- 
dient is the existence of magnetic perturbations which 
explicitly break local spin-rotational invariance (e.g., due 



to local spin-orbit coupling ||) and thereby induce sub- 
stantial hole pair binding. We then develop a mean-field 
theory of superconductivity based upon a phenomenol- 
ogy from our microscopic model. We emphasize that in 
our approach there are two, in principle different, en- 
ergy scales; one associated to the pairing of holes and 
another related to the phase coherence of the pairs (that 
establishes T c ). Basically the inhomogeneities induce a 
strong hole pairing, which in turn Josephson-tunnel co- 
herently between stripes, separated by insulating AF re- 
gions, phase-locking into a macroscopic supercurrent su- 
perfluid stiffness. Recently, a simple linear relation be- 
tween the superconducting transition temperature T c and 
the AF incommensuration S has been observed for the 
LSCO § and YBCO (§ high-T c compounds: k B T c oc 5, 
where ks is the Boltzman constant. We find that for this 
relation to hold we need a power law Josephson tunnel- 
ing. 

Microscopic Model. Our model Hamiltonian describ- 
ing the low energy dynamics of Cu02 planes is H = 
Ht-j + flinh , where the background Hamiltonian H t -j is 
the standard t-J model 

H t -J = -t ^ 4a C jc + J J2( Si ' S i ~ 4"'%) • 

(i,j),<r (id) " (1) 

For the inhomogeneous component, we take 



inh 



2^ sj, s' a s' p 



5J± 



with SJ± ^ SJ Z , representing the magnetic perturba- 
tion of a static local Ising anisotropy, locally lower- 
ing spin symmetry. Here are near-neighbor sites, 
while (a, 0) are two near-neighbor sites characterizing 
the bonds that are perturbed and where SU(2) spin- 
rotational invariance is explicitly broken. The network 
of perturbed bonds form mesoscopic patterns determined 
by the distribution of stripe segments. The spin-^ oper- 
ator Si = \&i C! T (J(J ,c ia , , the electron occupation number 
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hi = ct^c^ + c li c iii an d c la( c ia) creates (annihilates) an 
electron of spin a in a Wannier orbital centered at site i; 
r are the 3 Pauli matrices. This is a three-state model 
with the hopping constrained to the subspace with no 
doubly occupied sites. In the following, all energies will 
be measured in units of J. 

Our modeling strategy consists in assuming the exis- 
tence of an inhomogeneous mesoscopic skeleton of stripe 
segments and then exploring its consequences, mainly 
the competition between magnetism and superconduc- 
tivity. We do not address here the important problem 
of the formation and stability of this skeleton morphol- 
ogy. The origin(s) of "stripe segment" formation in high 
temperature superconductors is as yet unclear and sev- 
eral physical mechanisms could act cooperatively and be 
responsible for the generation of multiple length scales, 
among them: spin-orbit coupling, local Jahn- Teller dis- 
tortions induced by the hole, effective interactions com- 
ing from a multi-band Hubbard Model (HM) (including 
explicitly the oxygen and copper bands), oxygen buck- 
ling at the stripe, and other local magnetoelastic effects 
. Competitions between attractive short range forces 
and repulsive long range ones can certainly spontaneously 
break translational and/or rotational invariance in the 
Cu02 planes plE2|, but this is not necessarily the only 
mechanism. However, we show below that the mere exis- 
tence of appropriate local magnetic anisotropies is crucial 
for pair-formation. 

We start by showing that, as far as we could numeri- 
cally determine, only by including a local Ising perturba- 
tion such as -f/inh in Eq. [I] can a strong pairing of holes 
be obtained (see Fig. [I]). All the calculations were made 
using exact diagonalization in small clusters with peri- 
odic boundary conditions in all spatial directions. We 
studied one-dimensional (ID) chains up to 16 sites and 
8x2 clusters. Hereafter, we will view our clusters as 
simulating systems where the longer direction is perpen- 
dicular to the stripes. We investigated the system size 
scaling for the binding energy of two holes defined as 
Eb = (Wholes — -Eohoic) — 2(Ei hole — -Echoic) for sev- 
eral models in ID and 2D. Although for small enough 
systems the binding energies could be very large, they 
all seem to extrapolate towards no (or extremely small) 
binding in the thermodynamic limit, with the clear ex- 
ception of the inhomogeneous t-JJ z case. We have also 
studied several one-band HMs, but we could again not 
find definite binding. The t-J J z model, the only one un- 
ambiguously giving binding in the thermodynamic limit, 
is obtained by breaking spin-rotation symmetry in d near- 
neighbor bonds (a, (3), repeated with period P, by an 
amount 6J Z = 0, 6J±_ < in Eq. |l|. This t-JJ z model 
is a most natural way to induce a spin-gap. We have 
checked that the spin-gap is present for our t-JJ z model 
p3[ . The inhomogeneities forming the superstructure, 
which we impose by hand in the Hamiltonian, we term 
stripes. In Fig. we show the hole correlation function 



(g\no.rii\g), where \g) is the ground state of the system 
{{g\g) = 1). This correlation function gives information 
about the structure of the pair. It can be seen that as 
the hopping strength t is increased beyond a character- 
istic value the second hole jumps from one stripe to the 
neighboring one, starting from an initial configuration 
where both holes are in the same stripe for small t. This 
can be understood as a result of a length(time)-scale com- 
petition: the pair size exceeds the stripe width. 

To explore the nature of the binding, we have examined 
the canonical transformation of a t-J model from a one- 
band HM and traced what kind of perturbations would 
produce a t-J J z term. This corresponds to a term like 
—VJ2a fH,<T^j,a (hi first neighbors) in an extended HM, 
which may in turn arise, for example, from local mag- 
netoelastic (e.g. oxygen (un)buckling) jl4| or spin-orbit 
couplings. Note, again, that here this is a perturbation 
only at the stripes. This kind of anisotropy manifests it- 
self in two different ways in the t-JJ z model, enhancing 
both the —jhihj and the easy axis (Ising) terms of the 
Hamiltonian. The first one is an explicit pairing term for 
electrons. To see the relative importance of each term we 
have calculated the binding energy of a Hamiltonian like 
(|l|) but excluding the —jhihj term and including bonds 
with broken spin-rotational symmetry. This model corre- 
sponds to holes (with no spin) propagating in an antifer- 
romagnet, but not derived from a canonical transforma- 
tion of a one-band HM. We find that it still has binding, 
as should be expected. Thus, the easy axis exchange term 
is partially responsible for the binding energy. 

In order to understand this exchanged-based pairing 
mechanism, it is useful to explore some limiting cases. 
When the magnetic energy scales are the most relevant 
ones: ( J z = J + 5J Z , J± = J + SJ±) 3> t, it is easy 
to realize that, depending upon J Zl J± being smaller or 
larger than J, the holes will prefer to be in the stripes 
or between stripes (with no binding), respectively. The 
opposite limit, i.e. purely kinetic energy, leads to delocal- 
ized holes and no binding. The situation where J z < J 
and t is relevant corresponds to the intermediate regime 
where pairing is observed. Notice that pairing of holes 
does not necessarily imply that holes should share the 
same stripe, they can occupy neighboring ones (see Fig. 
^, upper panel), thus avoiding phase separation. Details 
of the charge confinement and pairing potentials from the 
(dynamic) spin-field profiles in the superlattice skeletons 
will be given elsewhere. 

Having demonstrated a minimal model for hole bind- 
ing, we have computed spin correlation functions in clus- 
ters of size N x x N y = N (N x = 8,N y = 2). Here, we 
simulate the stripes by including an anisotropic 6J± < 
in one y-bond with P — 4; the rest of the bonds, including 
all the x-bonds, were not changed from the background 
t-J model (see inset, Fig. ^). We cannot perform scaling 
on this size of inhomogeneous system, but the binding 
energy is still considerable. We have included up to 6 
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holes. In the case of four holes (the one more relevant to 
the stripes in the underdoped regime for cuprate oxides) 
and small t (< J) the holes bind in pairs on each site of 
the inhomogeneous bond (see Fig. ||, lower panel). 

In Fig. H we show the spin-structure factor 
function S(k = {k x ,k y )) defined as: (SkS-k) = 
(1/N) Y^,i j exp(ik.ri) (g|Sj.Sj + i|ff). This function corre- 
sponds to the observable in the elastic neutron scatter- 
ing experiments. For t small, only one peak occurs with 
k y = 7T, corresponding to two essentially uncorrelated 
AF domains, isolated from each other by the pinned hole 
wall. As t increases (t/J> 2, near the accepted set of val- 
ues of the 2D t-J model for cuprates fTsfl ), the holes gain 
kinetic energy by visiting the first neighbor sites around 
the anisotropy region, but still bind together. The effec- 
tive width of the pair thus increases to two sites. Mag- 
netic energy is then gained if the two domains shift their 
staggered magnetization by n. We suggest that these 
0(t 2 ) processes are responsible for the incommensuration 
(8) in SQs.) observed in the experiments 0. This 8 is the 
inverse of twice the period P of the stripes. In this pic- 
ture the incommensuration is a consequence of the holes 
and their kinetic energy, and is a property of the ground 
state. Basically, it results from the competition between 
hole delocalization and magnetic fluctuations. This con- 
trasts with some other proposed explanations, where 8 is 
a magnetic thermodynamic property fll6|| . It is interest- 
ing to note that this incommensuration arises even in the 
homogeneous t-J model although for different values of 
t. This suggests that the experimentally observed mag- 
netic properties are already present in a homogeneous t-J 
model, but in order to obtain binding of holes appropiate 
inhomogeneous terms must be included. 

When more than four holes are added to the system, 
but only two bonds are perturbed, S(k) changes qualita- 
tively. Instead of showing an incommensurability around 
k = Q = (71-, 7r), it has a broad peak at k = (0, w). In this 
case the extra holes are delocalized in the middle of the 
AF space between stripes. This suggests that when the 
stripes reach their minimum separation, extra holes are 
responsible for the experimental increase and ultimate 
disappearance of the incommensurability. 

Model of Josephson Spaghetti. It is important to relate 
the above discussion to the experimental evidence for the 
incommensurate neutron scattering peak, seen in LSCO 
(e.g., ||) and YBCO compounds 0. In both of these 
cases a simple linear relation between T c and the peak 
incommensuration 8 near Q (or peak width in YBCO) is 
obeyed f||§. Namely, 

k B T c = hv*S . (2) 

The anomalously low velocity values for v* depend on 
the compound Q. These velocities are independent of 
the carrier concentration and the only doping (x) depen- 
dence entering Eq. is through 8(x). 



An interpretation of this relation is to connect pos- 
sible superconductivity mechanisms to the existence of 
the fluctuating stripes. Here we focus on the simple 
proportionality between T c (x) and a doping dependent 
length £{x), determined from the neutron scattering: 
T c (x) oc l/£(x), £(x) = l/8(x). We consider how the 
Josephson tunneling of pairs between stripe segments can 
produce the relation between the phase ordering transi- 
tion temperature T c and the typical length £(x). The 
stripe-stripe distance r is a random quantity due to in- 
trinsic mechanisms as well as disorder and/or crystal 
imperfections Therefore, we will assume that 

the mean-field transition temperature depends upon the 
Josephson coupling (J(r)), averaged with some probabil- 
ity distribution of stripe separations. 

Our model Hamiltonian of random stripe separation 
and associated inter- and intra-stripe random Josephson 
coupling (see Fig. |J) is 

H = ^2 Jij exp[i (</>,; - <f>j)} , 

ij 

J ij =J(r ij )=t Q /rf j , (3) 

where the summation is taken over the coarse-grained 
regions i — 1, ••• ,JV with well-defined phases, labeled 
4>i = 4>{ri) and becomes zero eventually at large dis- 
tances. Next, we will assume some probability distribu- 
tion P{r) for the stripe-stripe distance. For simplicity 
we will take the "box" distribution P{r) centered around 
1=1/8 and with finite width a = v£, where v = 0(1) is 
a parameter. P(r) = C, for £ — a < r < £ + a, and zero 
otherwise. Here we have simplified to one length scale 
for both a and £. The normalization constant in 2D is 
C = [4tt£o\~ 1 . In this model one easily finds 

<J(r)> = f <PrP(r)J(r) = ^^a^ , 
J 2 - a 

(r) = ^a 2 e , (4) 

where the constants a\, a-i are 0(V). Thus, for a = 1, we 
obtain the experimentally observed relation 

T c (x)^(J(r))^[(r)}-'=S(x) . (5) 

We have examined a variety of distributions P(r) and 
functional dependences for Jy; Eq. || with a = 1 is 
the only one reproducing the experimental data (at our 
mean-field random Josephson coupling level. Implicit in 
J(r) is the exponential cutoff at lengths much larger than 
the stripe-stripe distance. This cutoff is necessary to have 
a well defined thermodynamic limit but is not important 
for short length scales). The screening mechanism (mag- 
netic, elastic fluctuations, etc.) responsible for this form 
requires detailed microscopic modeling [ p"7| . The present 
model does not allow us to determine the magnitude of 
v* without making specific assumptions about parame- 
ters such as to. 
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In conclusion, we have presented a microscopic model 
that captures the essential magnetic and pairing proper- 
ties of high-temperature cuprate superconductors. Pair- 
ing of holes is a consequence of the existence of an AF 
background. (Analogous scenarios in other broken sym- 
metry backgrounds, e.g. doped charge-density-wave bis- 
muthates, are likely.) Crucially, however, the glue is pro- 
vided by magnetic inhomogeneities whose precise origin 
remains to be unraveled, although it seems fundamen- 
tal that these perturbations should locally break spin- 
rotational invariance. This pairing mechanism is kinetic 
exchange-interaction based and involves a competition 
between Ising and XY symmetries. We emphasize that 
the pair-binding occurs only for intermediate strengths of 
t and (local) J z . We also introduced a phenomenological 
model and scenario for the macroscopic superconductiv- 
ity based upon coherent Josephson-tunneling of pairs of 
holes between these magnetic inhomogeneities in a meso- 
scopic liquid-crystal-like Jll|] skeleton. We have shown 
that this approach is able to recover the magnetic incom- 
mensuration 8 and its experimentally observed relation 
to T c (x). Finally, we note that we have assumed static 
magnetic inhomogeneities. The case where the broken 
spin-symmetry follows the hole is also interesting. Else- 
where, we will discuss this generalization of coupling the 
inhomogeneity selfconsistently to dynamic holes. 

Work at Los Alamos is sponsored by the US DOE un- 
der contract W-7405-ENG-36. 
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FIG. 1. Schematic representation of the superconducting 
ground state. AF zones between stripe segments are col- 
ored in blue. Red spins near the stripes represent easy-axis 
Ising-like links. Gray ellipses characterizes the bound pairs 
(holes in red). Note the zig-zag alignment of the holes, and 
the AF domains 7r-shifted at each side of the stripes. 



FIG. 2. Correlation function for 2 holes in a 16x1 chain 
with d — 3, P — 8, J = 1 and SJ± — —0.9 (top); and for 
4 holes in 8x2 with d = 1, P = 4 and the other parameters 
as before (bottom). Lines in the insets show the perturbed 
bonds. In 16x1, as t grows the pair switches from occupying 
one single stripe, to two neighboring ones. In the (2 stripes) 
8x2 case, the holes remain at both ends of the anisotropic 
bond, but as t is increased they form an effective two-site 
stripe because of transverse spin fluctuations. The magnetic 
energy is lowered by 7r-shifting the AF domains on either side 
of the stripe. 



FIG. 3. Spin structure factor for a t-JJ z ladder (8x2) with 
two SJ± = —0.9 bonds in the {/-direction (see inset). The 
incommensurability appears only for t larger than a critical 
value. When 6 holes are added to the system the double peak 
disappears and is replaced by a broad one around k = (0, 7r). 



FIG. 4. Schematic Josephson coupling between an as- 
sumed distribution of stripe segments. For the incommen- 
suration 8 to be observed along crystallographic (1,0) and 
(0,1) directions, the stripe-stripe distances must have average 
(r) ~ I = 1/8. (J) will be determined by the probability 
distribution P(r) that determines the statistics of inter-stripe 
distances. Physically it is clear that P(r) should be centered 
near £, with some width arising from the meandering of stripes 
(see text). 
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